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ABSTRACT 

We determine the entire electrical current density vector in a geometrical 3D volume of the in- 
ner penumbra of a sunspot from an inversion of spectropolarimetric data obtained with Hinode/SP. 
Significant currents are seen to wrap around the hotter, more elevated regions with lower and more 
horizontal magnetic field that harbor strong upflows and radial outflows (the intraspines) . The hori- 
zontal component of the current density vector is 3-4 times larger than the vertical; nearly all previous 
studies only obtain the vertical component J z and thus strongly underestimate the current density. 
The current density J and the magnetic field B form an angle of about 20° . The plasma /3 at the km 
level is larger than 1 in the intraspines and is one order of magnitude lower in the background com- 
ponent of the penumbra (spines). At the 200km level, the plasma /3 is below 0.3 nearly everywhere. 
The plasma /? surface as well as the surface optical depth unity are very corrugated. At the borders of 
intraspines and inside, B is not force-free at deeper layers and nearly force free at the top layers. The 
magnetic field of the spines is close to being potential everywhere. The dissipated ohmic energy is five 
orders of magnitudes smaller than the solar energy flux and thus negligible for the energy balance of 
the penumbra. 

Subject headings: Methods: numerical, observational - Sun: magnetic topology, sunspots - Techniques: 
polarimetric 



1. INTRODUCTION 

The study of the stability or dynamics of penum- 
bral filaments requires the knowledge of the Lorentz 
force, i.e. an accurate determination of the electri- 
cal current density vector J and the magnetic field 
vector B. The estimation of the energy budget dis- 
sipated by electrical currents obviously requires the 
determination of J. However, it is not trivial to 
reliably derive the electrical currents from observa- 
tional data. Previous attempts were aimed almost 
exclusively at the determination of only the vertical 
component J z with different degree of sophistication in 
the analysis of the data. The bulk of such studies was 
based o n magnetogram s dDeloach et al.lll984t lHagvard 
1988[ iHofmann etHI 119881 119891: IC anfield et al 



1993; 



1992; Ide la Beauiardiere et al] 19931: iLeka et al 



Metcalf et a 
iWang et al 



1994; 

Gary fe DemoulirJ 



1994, 



Ivan Driel-Gesztelvi et al 



19941: iZhang fc Wand 119941: 
19951: ILi et al.1 119971: IGao et al.1 
20081) . More accurate estimates of J z stem from 
state-of-art spectropolarimetric observations and the 
application of inversion techniques. In case of Milne- 
Eddington (ME) inversions (e.g. Skuman ich fc Lite's! 

119871 lLagg et all 120041 ) . B is obtained at an average 
optical depth. Among the recent wor ks calculating 



under ME approximation one finds IShimizu et al 



2009), iVenkatakrishnan fc Tiwaril (|2009h . and ILi et al 
(2009). A n alternat iv e dete r mination of J 7l is ob - 
taine d bv iBalthasarl (1200611 Uurcak et all (|2006l) . 
and iBalthasar fc Gomorvl (|2008l ) employing SIR 
(Stokes Inversion based on Re sponse functions, 
iRuiz Cobo fe del Toro Iniestal |l992), which delivers 
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the stratification of B in an optical depth scale. 

There have been several attempts to obtain Jhor, 
the modulus of the horizontal component of J, impos- 
ing app r oxima tion s to the magnetic field distr ibution: 
Ui et al.l (j2003| ) and lGeorgoulis fe LaBontel (|2004f ) obtain 
a lower limit of Jw assum ing a field- free configuration; 
iPevtsov fc Peregudl (|1990l ) derive Jhor imposing cylindri- 
cal symmetry of B in a sunspot. Without making any 
hypothesis on B, the determination of the three com- 
ponents of J requires the knowledge of the entire mag- 
netic field vector in a geometrical 3D volume, i.e. the 
determination of a geometrical height scale is manda- 
tory. While commonly in quiet Sun the transforma- 
tion from an optical depth scale to geometrical heights 
is done by assuming h ydrostatic equilibrium (see, e.g. 
iPuschmann et al"1l2005[ ). this is not justified in the mag- 
netized penumbra. In this case the force balance must 
include magnetic forces which requires to calculate the 
horizontal and vertical spatial derivatives of B. A first at- 
te mpt to empirica ll y deri ve 3D vector currents was done 
bv iSocas-Navarrol (|2005| ). w ho determined a geometri- 
cal height scale (following iSanchez Almeida! 120051 ) by 
imposing equal total pressure between adjacent pixels, 
although neglecting the magnetic tension in the Lorentz 
force. 

For an accurate determination of J in the present pa- 
per, we take advantage of the 3D geometrical model of a 
section of the inn e r pen umbra of a sunspot described in 
IPuschmann et all (|2010D [hereafter, Paper I]. We use ob- 
servations of the active region AR 10953 near solar disk 
center obtained on 1 st of May 2007 with the Hinode/SP. 
The inner, centerside, penumbral area under study was 
located at an heliocentric angle 6> = 4.63°. To derive the 



2 



Puschmann et al. 




12 3 4 

x [Mm] 



FlG. 1. — Upper panel: B z (in color) and B^ OI (arrows) at 
z = 200 km. The size of the arrows ranges from 1220 to 1900 G. 
Contour lines enclose areas with J> 120 mA/m 2 . Lower panel: J z 
(in color) and J^ OI (arrows) at z = 200 km. The size of the arrows 
ranges from 1 to 365 mA/m 2 . Contour lines correspond to B z equal 
to -650 G (solid lines) and to -450 G (dashed lines). 

physical parameters of the solar atmosphere as a func- 
tion of continuum optical depth, the SIR inversion code 
was applied on the data set. The 3D geometrical model 
was derived by means of a genetic algorithm that mini- 
mized the divergence of the magnetic held vector and the 
deviations from static equilibrium considering pressure 
gradients, gravity and the Lorentz force. For a detailed 
description we refer to Paper I. 

2. RESULTS AND DISCUSSION 

The current density vector J=(V x §)/hq has been 
calculated for each pixel in the FOV analyzed in Paper 
I at geometrical height layers between and 200 km, 
i.e. in a volume of (4.2 Mm x 5.6 Mm x 0.2 Mm) in the 
inner penumbra of a sunspot placed close to the disk 



TABLE 1 

Mean values of the moduli of the vertical and horizontal 
component of j as well as the mean value of j at three 
different height layers. 



z 

[km] 


<| Jz |> 

[mA/m 2 ] 


< Jhor > 

[mA/m 2 ] 


< J > 
[mA/m 2 ] 





30 ±20 


129 ±50 


135 ±59 


100 


26 ± 9 


103 ± 36 


109 ±39 


200 


25 ± 8 


88 ±33 


93 ±35 
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center. We decompose B = B z e z + £?hor where e z de- 
note the unit vector in vertical direction. In a similar 
way, we write J = J z e m + Jhor- Through the paper B, 
-Bhor, J and Jh or will denote the modulus of the cor- 
responding vectors. The upper panel of Fig. [T] shows 
B z (colored background) and -Bhm (black arrows) at the 
top layer (200 km). Areas in red color correspond to 
regions with lower and more horizontal magnetic field 
(intraspines. lLites et al.lll993l ) which are hotter, elevated 
and harbor strong upflows and radial outflows (see Pa- 
per I). These regions can be interpreted as the embedded 
nearly horizontal flux tubes of the uncombed scenario 
(see e. g. lSolanki fe Montavodll993tlSchlichenmaier et all 
Il998alfbl : Martinez Pilletl 120001: IBorrero fc Solankil 120101 
and references therein). We will call spines to regions 
with a more vertical and intense magnetic field (they cor- 
respond to the background component in the uncombed 
scenario). White contour lines enclose areas with signifi- 
cant J values (larger than 120 mA/m 2 ), located predom- 
inantly along the borders of the intraspines. The lower 
panel of Fig. [T] shows J z (colored background) and Jhor 
(black arrows). White contour lines correspond to B z 
values equal to -450 (solid) and -650 G (dashed), respec- 
tively. The electrical currents wrap around the nearly 
horizontal flux tubes (intraspines): For the majority of 
the intraspines, J z shows positive values at the upper 
(larger Y-coordinate) borders of the filaments and nega- 
tive values at the lower borders. The orientation of the 
black arrows indicates that the currents circumvent the 
flux tubes in most cases. 

In the left panel of Fig. [2] we plot the histogram of 
J z evaluated at the top layer (200 km) in black. To 
check the significance of the results, we evaluate a sim- 
ulated electrical current density distribution J orr . The 
random vector field J err deviates from zero only by a 
Gaussian noise with a a equal to the estimated error at 

each pixel for each component of J; a was obtained from 
an error propagation of the uncertainties of B. The his- 
togram of J crr , z is represented in red. Frequently, the 
z component of the current density is evaluated from 
the results of ME inversions; to determine the reliabil- 
ity of such results we performed a ME analysis of our 
data set, and we calculated the vertical current density 
Jz,me- The histogram of J Zj me is plotted in blue in the 
left panel of Fig. [2J Note that the J z calculated from 
our 3D geometrical model is significantly larger than the 
Jz.me, however both are clearly above the calculated 
uncertainties. The ME inversion delivers the magnetic 
field just at an average optical depth a round logr=- 
1.5 (|Ruiz Cobo fc del Toro Iniestal Il99l) . As we have 
seen before, larger electrical currents appear at the bor- 
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Fig. 2. — Left panel: Histogram of J z at 200 km (black) and histogram of J z retrieved from a Milne-Eddington inversion (blue). Middle 
and right panel: histograms of Jh or an< i J evaluated at the same height layer. Red lines represent the corresponding error distributions. 

ders of structures with lower magnetic field and smaller 
optical depth. Consequently, the ME inversion probes 
higher layers above these structures resulting in a smaller 
derivative of the magnetic field as the one obtained from 
the 3D geometrical model. In the middle and right panel 
of Fig. [3 we present the histograms of Jh or and J, respec- 
tively. In red we show again the estimated uncertainties. 
As mentioned before, in most observational studies of 
electrical currents in solar active regions only the verti- 
cal component has been measured. We find Jhor to be 
about four times larger than J z . Note that the horizon- 
tal component is clearly predominant, with a distribu- 
tion practically equal to the one of J. Works estimat- 
ing Jhor by making certain h ypotheses on the B config - 
ura tion reach similar result s: Peytsov fc Peregudl (|1990f > 
and lGeorgoulis fc LaBontd (120041) fo und Jw ~ 2-3 times 
larger than J, . while Ui et al.l (j2003l ) report on Jh or val- 
ues of about one or two orders larger than J z . All studies 

estimating only J z strongly underestimate J present in 
sunspot penumbrae. 

Table 1 summarizes the values of < J z |>, < Jh or >, 
and < J > at three different heights. The errors given in 
the table have been calculated by an error propagation 
from the distribution of J crr described above. The electri- 
cal current density decreases with height and at all layers 
the horizontal component is on average about four times 
larger than the vertical one. While < Jh or > decreases 
by 32% between the and 200 km level, <| J z |> shows a 
small decrease with height of about 17%. The large un- 
certainties at the km level stem from the corresponding 
uncertainties of the determination of the magnetic field 
due to the low sensitivity of the visible lines used in this 
study. 

The upper left panel of Fig. [3] shows the histogram 
of the angle ip formed by B and J at the Okm height 
layer (black). In a force-free configuration, B and J 
are parallel to each other. In our volume ip is of about 
20°. Since the uncertainty of ip is large for pixels with 
small J, we show also the histogram of ip for pixels with 
J(z = 200) > 120 mA/m 2 (red line). The difference in 

the orientation of B and J is mainly due to a difference 
in their inclination from the vertical 7 and not due to 
a difference in their azimuth <f> as the upper right and 
lower left panel of Fig. [3] indicate: Both vector fields are 
basically in the same vertical plane. Consequently, the 
horizontal component of the Lorentz force must be larger 
than the vertical one. In fact < ( J x B)hor > amounts 
to 1.03 millidynes/cm 3 , whereas < (J x B) z > shows 



Fig. 3. — Upper left panel: Histogram of the angle ip between J 
and B at z = km. Upper right panel: Histogram of the difference 
between the inclination from the vertical of B and J at z = km. 
Lower left panel: The same for the difference in azimuth. Lower 
right panel: Energy dissipated by electrical currents in units of 
10 solar flux (between and 200 km, black; between -225 and 
200km, blue). In all panels the red line corresponds to angles and 
energies evaluated for those pixels with J(z = 200) > 120 mA/m 2 . 

a value of 0.61 millidynes/cm 3 only, both evaluated at 
z = 0km for pixels with J(z = 200) > 120 mA /m 2 . This 

m 



is in agreement with IPevtsov fc Peregudl (|1990f) who also 
found the horizontal component of the Lorentz force be- 
ing larger than the vertical one. 

The lower right panel of Fig. [3] shows the histograms 
of the energy dissipated by the electrical currents (ohmic 
energy) in the region under study integrated between the 
and 200km level (black). The ohmic energy has been 
calculate d using the longitud i nal el ectrical conductivity 
following IKopeckv fc Kukfinl (|1969D and is negligible in 
the inner sunspot penumbra compared with the solar 
flux. The energy is dissipated only at the borders of 
the horizontal tubes, where J reaches significant values 
(red line). In Paper I we found that the main contri- 
bution to the energy flux carried by the ascending mass 
(convective energy) stems from layers below -75 km. The 
resulting convective energy flux, integrated from -225 to 
200 km, reaches values of up to 78% of the solar flux, and 
thus would be sufficient to explain the observed penum- 
bral brightness. However, this result has to be taken with 
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Fig. 4. — Map of the plasma f3 at z = 0km. Contour lines cor- 
respond to B z equal to -650 G (solid lines) and to -450 G (dashed 
lines) . 

care since the physical parameters at layers below -75 km 
entering this calculation result from extrapolations. If all 
the layers between -225 and 200 km are included in the 
calculation of the ohmic energy, the resulting values (well 
below 10 -4 Fq) are still negligible for the energy balance 
(blue line in Fig. 

The plasma f3, defined as the ratio between the gas 
pressure and the magnetic pressure, is an important pa- 
rameter to study the relevance of the Lorentz force in 
the dynamics of sunspot penumbrae. Fig. U shows the 
plasma /3 at a height of km. In our inner penumbral 
area, the /3 = 1 surface is strongly corrugated. Note that 
in the intraspines (areas enclosed by white contour lines) 
the plasma j3 is clearly larger than 1, being one order of 
magnitude lower in the spines. At the top layer (200 km) 
not presented here, the aspect is similar, although all 
values have decreased by a factor ~ 5. 

3. CONCLUSIONS 

The application of a genetic algorithm on the optical 
depth model retrieved from a SIR inversion of spectropo- 
larimetric Hinode/SP data allowed us to construct a 3D 
geometrical model of a section of the inner penumbra 
of a sunspot (see Paper I). The resulting model has, by 



construction, a minimal divergence of B and a minimal 
deviation of the static equilibrium. This allows us for 
the first time the accurate determination of the entire 
3D electrical current density vector. At the 0'.'32 reso- 
lution of the Hinode/SP, the electrical current density 
appears to be spatially resolved in the penumbra: Sig- 
nificant J values (larger than 120 mA/m 2 ) are located 
at the borders of the intraspines. The electrical currents 
seem to wrap around the nearly horizontal flux tubes. 

We find an horizontal component of J about four times 
larger than the vertical one. This is in agreement with 
earlier works estimating J^or by assuming simplifying hy- 
potheses on the B distribution. All works 1 only evaluat- 
ing J z clearly underestimate J. 

The magnetic field at lower layers is not force-free at 
the borders^of the intraspines: There, the angle formed 
by B and J is about 20°. The difference in the orienta- 
tion of B and J is mainly due to a different inclination 
with respect to the vertical of both vectors, being nearly 
negligible the difference in azimuth, leading to a domi- 
nant horizontal component of the Lorentz force, directed 
towards the central axis of the intraspines, something 
that helps them to maintain the internal force balance. 
At the km level, the plasma (3 is strongly corrugated, 
being larger than 1 at the borders and inside the in- 
traspines and one order of magnitude lower inside the 
spines. In the latter, the Lorentz force can hardly be 
balanced by the pressure gradient or weight of the mate- 
rial: consequently the magnetic field configuration of the 
spines in the inner penumbra must be close to a force- 
free configuration. Furthermore, as in these areas J is 
relatively small and B is large, the field must be closer 
to being potential. Following the same reasoning, our 
results show that at the highest layers the plasma j3 is 
nearly everywhere below 0.3 and thus the field must be 
closer to a force- free configuration almost everywhere and 
to a potential one in the spines. 

The dissipated ohmic energy is clearly negligible being 
5 orders of magnitudes smaller than the solar flux. 
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